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We explore the properties of weakly bound bosonic states in the strongly interacting regime. Combining 
a correlated-Gaussian (CG) basis set expansion with a complex scaling method, we extract the energies and 
structural properties of bosonic cluster states with N < 6 for different two-body potentials. The identification 
of five- and six-body resonances attached to an excited Efimov trimer provides strong support to the premise of 
Efimov universality in bosonic systems. Our study also reveals a rich structure of bosonic cluster states. Besides 
the lowest cluster states which behave as bosonic droplets, we identify cluster states weakly bound to one or two 
atoms forming effective cluster-atom "dimers" and cluster-atom-atom "trimers." The experimental signatures 
of these cluster states are discussed. 



The understanding of the universal nature of low-energy 
few-body body physics is a fundamental prerequisite for the 
development of effective many-body descriptions in the ultra- 
cold regime. In two-component Fermi gases, the character- 
ization of two-body physics in terms of a single interaction 
parameter (the scattering length) is at the heart of our under- 
standing of the BCS-BEC crossover. In bosonic systems, the 
underlying few-body physics is enriched by Efimov physics 
which leads to the formation of a series of trimers which ac- 
quire peculiar properties such as a borromean nature and a 
discrete scale invariance UEJ. Recent experimental devel- 
opments have opened a new era in the exploration of Efimov 
phenomena in ultracold gases [3-6 1. While this richness and 
complexity of three-boson systems have led to a new level of 
understanding of universal few-body physics, they pose im- 
portant questions as to whether the low-energy behavior of 
larger bosonic systems can be understood and characterized 
within a simple universal framework. 

The natural starting point for addressing this question is the 
exploration of Efimov and universal phenomena in increas- 
ingly larger systems. For N > 3, a quantitative analysis of 
universality is significantly more challenging because it usu- 
ally requires the study of resonances rather than bound states, 
since bound states occur in a regime where the applicabil- 
ity of universality is questionable. However, in the last few 
years, tremendous progress has been achieved in the under- 
standing of universality at the four-boson level flT HTOl . In par- 
ticular, a universal regime has been recognized in which the 
four-body physics is determined by the underlying two- and 
three-body physics. Recently, some of these early four-body 
predictions have been quantitatively verified with remarkable 
accuracy ifTTl . Despite these important advances, there still 
remain controversies regarding the scope of universality, the 
role of the four-body parameter, and the nonuniversal cor- 
rections ifTZlfPfl . For N > 4 systems, the applicability of 
universal theory is even more debatable because of the lack 
of theoretical or experimental evidence of universal behav- 
ior. A natural continuation to the four-body predictions of 
Ref. E] would indicate that universality extends to larger 
clusters, whose behavior follows the same Efimov discrete 
scale invariance and is only controlled by two and three-body 
physics. Thus, for each Efimov state (N — 3), there would 
be a series of iV-body states (or resonances) with N > 3 as- 



sociated to it forming an Efimov family. A recent study |fT31 , 
based on such premise of universality, has characterized some 
of properties of the lowest weakly bound cluster states up to 
N < 13. While this study provides key predictions to be theo- 
retically and experimentally explored, it leaves open questions 
such as the validity of the universality hypothesis and the ex- 
istence of additional weakly bound cluster states. 

In this article, we address these questions through the analy- 
sis of the structure of the strongly interacting few-boson spec- 
tra. First, we explore the universal regime and we identify 
five- and six-boson resonances tied to an excited-Efimov state. 
These resonant states represent small bosonic droplets and are 
in good qualitative agreement with the prediction of Ref. lfT31 . 
Our results provide much-needed support to the premise of 
universality and discrete scale invariance in few-boson sys- 
tems. We extend our analysis to the exploration of different 
types of cluster states formed in the lowest Efimov family, 
which is strongly modified by nonuniversal corrections. For 
a range of model potentials, we find that cluster states (N — 
3, 4, 5) are likely to bind weakly to atoms forming effective 
cluster-atom "dimers" and cluster-atom-atom "trimers." We 
identify one of these cluster-atom-atom "trimers" as a reso- 
nance that appears energetically slightly below the lowest Efi- 
mov trimer and it can be qualitatively described as an Efi- 
mov trimer formed by an Efimov trimer and two atoms. This 
state is one of the simplest bizarre cluster structures mathe- 
matically proposed lfl6l : similar cluster-atom-atom structures 
are expected for larger systems. Our studies reveal an intricate 
structure of bosonic cluster states and provide estimates of the 
./V-body resonant positions relevant to experiments. 

Our starting point is the few-boson Hamiltonian, 

i i<j 

where m is the mass of the bosons and ry is the inter- 
particle distance between particles i and j. The two-body 
model potential takes form V(r) — Vb(exp[— r 2 /(2d 1 )] — 
a cxp[— 2r 2 /d^]), where do is the interaction range, and Vo 
and a are tuned to change the shape and the scattering length 
of the potential. We consider two qualitatively different po- 
tential: purely attractive (Vo < and a < 1) and attractive 
with repulsive core potentials (Vq < and a > 1). While 
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the studies considered a —20 < a < 20, most of the results 
will be presented for two cases: a purely attractive interaction 
V a (a — 0) and an attractive potential with a soft repulsive 
core V r (a = 2 and Vo < 0). The range r and the energy 
E sr = ti 2 / (mrg) are the typical energy and length scales that 
characterize the interactions. The universality regime is char- 
acterized by energies <C E sr and length scales (character- 
izing scattering length and clusters sizes) £ S> r^. 

To describe the bound states, we use a CG basis set ex- 
pansion ifTTl fTSll that has been very successful in describ- 
ing bosonic and fermionic systems with short-range inter- 
actions. In our implementation, the eigenstates of a sys- 
tem are expanded in the set of CG basis functions in which 
the center-of-mass coordinate has been removed, the rela- 
tive angular momentum is zero and the parity is positive 
(L 7 ' = + ). Each basis function is a symmetrized prod- 
uct of Gaussian functions, each of which depends on one 
of the N(N — l)/2 interparticle distances and can be writ- 
ten as ijjp{x) = exp (- £V ^fj x « ■ where x = 
{xi,X2, ...,xjv-i} is a set of (relative) Jacobi coordinates. 
The parameters A^j that characterize the Gaussian function 
widths are selected and optimized using a stochastical varia- 
tional method IfTTl . The convergence of the results is carefully 
analyzed by increasing and reoptimizing the basis set. Typi- 
cal calculations for N > 3 include up to 3500 optimized and 
fully symmetrized basis function. To explore the structure of 
the few-body states, we extract the pair-distribution function 
defined as 47rr 2 P/v(r) = (^N\^(rij — r)|^Ar). where r,j is 
the interparticle distance between particles i and j. 

To study resonances, we use the complex-scaling method 
(CSM) H2HID- In the CSM, all coordinates are rotated as 
r — > re 10 by a transformation U (9). The wave function of the 
resonance is square-integrable in these rotated coordinates and 
can be expanded in the same square integrable basis functions 
that describe bound states: 

* e (x) = ^)*(x)=5^C i (^ i (x) 1 (2) 

i 

where x = {xi, X2, xjy-i} is a set of (relative) Jacobi 
coordinates. The wave function ^e(x) is a solution of the ro- 
tated Hamiltonian He = U(8)HU(9)^ 1 with complex energy 
Eg = Eji — iT/2, where T is associated with the width of the 
resonance. 

Universal Droplets - Figure [T] summarizes the energies of 
the universal cluster states for N — 3 — 6. The results are 
presented as a function of the relevant universal parameters: 
1/na and E/E%, where k = \J mE^ / h 2 is the three-body 
parameter, and E% is the binding energy of the Efimov trimer 
at unitarity (a = oo). The five- and six-body results are ob- 
tained from an analysis of resonances attached to an excited 
Efimov family. These resonances are observed for both the V a 
and V r potentials. At unitarity, we obtain Eg f=s 10.1 £3 in 
close agreement of Ref. 031 predictions of E£ sa 10.4(2) E% . 
For N = 6, we find that E% w W.3E% which is slightly 
lower than the predictions from Ref. US) of E% sa 18.4(2)£^. 
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FIG. 1: (color online) Universal cluster energies as a function of 
inverse scattering length. The energies for N = 3 and N = 4 were 
obtained in previous studies 1 8|. The N = 5 solid curve and symbols 
correspond to predictions from an excited resonance for V r and V a 
respectively. The circles and squares are the N = 6 predictions for 
V r and V a respectively and the solid line is a guide to the eye. Dashed 
curves correspond to the TV = 5, 6 predictions from Ref. 1 15 1. Inset: 
Width of the five-body resonances for V r (red) and V a (blue). 



Interestingly, the energies of the universal states at unitar- 
ity scales roughly with the number of trimer configurations 
the cluster supports, i.e., 1, 4, 10, 20 for N = 3, 4, 5, and 
6. To further verify the agreement between different model 
potential-predictions, we analyze the pair distribution func- 
tion at unitarity (see Fig. |2j. The good agreement between 
three different predictions of P5 and P§ demonstrate the uni- 
versality of such few-body resonances. 
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FIG. 2: (color online) Pair distribution function of the universal five- 
body state. Here kn = \JmE^/h 2 , where E^, is the binding 
energy of the iV-body droplet. Solid and long-dashed curves are 
the prediction from the V a and V r , potential, while the short-dashed 
curve is the prediction from Ref. [15|. Inset: pair distribution func- 
tion for the six-body state. Same curve style convention as the main 
figure. 

The widths of the four- and five-body resonances depend 
strongly on the open decay channels. In our calculations, res- 
onances of the first excited Efimov family can only decay into 
the lowest Efimov family, i.e. into decay channels that are not 
quantitative in the universal regime. Therefore we expect the 
width of the resonance to be more sensible to non-universal 
corrections. For example, our analysis of the first-excited 
four-body resonance leads to a T^ b w O.IE^ , which is a fac- 
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tor of 3 larger than the predicted-converged width Ref. ifTTI 
(similar deviations were reported |QT|). The width of the first 
excited five-body resonance is presented in Fig.[TJfor both V a 
and V r . The clear differences between the different Y pre- 
dictions illustrate the importance of the nonuniversal correc- 
tions for resonances belonging to the lowest-excited Efimov 
families. However, as in the four-body case, we expect that 
the width extracted from the first-excited five-body resonance 
provides a correct order of magnitude estimate of the widths 
in the universal limit. We also estimate ~ O.SEg . 

The description of these five- and six-body resonances is 
extremely challenging since it entails an exploration of ener- 
gies that are a factor of e 27r / s ° ss 515 smaller than the en- 
ergies of the lowest fragmentation threshold. To obtain an 
accurate representation of these resonances, we carry out a 
numerical procedure inspired by previous implementations of 
the SVM+CSM |20]|2T|. First, we generate a basis optimized 
to describe a wide range of energies. We introduce an ex- 
ternal trapping potential V ex t(r) = mw 2 r 2 /2 in this part of 
the optimization to reduce the number of states described at 
each fragmentation threshold. As the basis is increased, the 
trapping frequency is reduced and finally set to zero to ex- 
tract bound states and resonances. The energy trajectories as 
a function of the basis set size show clear plateaus around the 
energy resonance that helps us identify it providing first esti- 
mates of the real part of the energy. To improve the descrip- 
tion of the resonances, a second optimization procedure is car- 
ried out involving only states in an energy window around 
the resonance. Once the optimization is complete, we cal- 
culate the spectrum as a function of 9 to identify the energy 
and the width of the resonance. The large number of avoided 
crossings between the resonant state and other states makes 
it particularly challenging to quantitatively estimate the six- 
body energy (see symbols in Fig. [TJ. Taking into account the 
numerical uncertainties and the potential nonuniversal correc- 
tions in the first excited Efimov family, we estimate a 10% 
uncertainty in the energies of resonances and the positions of 
the resonances discussed below. 

Excited cluster states and nonuniversal corrections.- Next, 
we analyze the formation of bound states for 3 < N < 6 in 
the lowest Efimov family. This study illustrates the structure 
of the spectrum in the strongly interacting regime by identify- 
ing other types of bound states. It also addresses the issue of 
nonuniversal corrections in the regime \a\ > ro that are par- 
ticularly important for understanding 133 Cs and 7 Li Efimov 
experiments in which five- and six-body resonant phenomena 
are expected to occur at \a\ < 3r v dw> where r v d w is the Van 
der Waals length. The general structure of the bosonic spec- 
trum is shown in Fig. [3] These results correspond to the po- 
tential V a . This structure, although it changes quantitatively, 
remains qualitatively the same for a range of model poten- 
tials. The lowest TV -body state is analogous to the universal 
states of Fig. [T] However, the energy of the lowest TV-body 
states grows very fast with the number of particles, implying 
that nonuniversal corrections increase with N (in agreement 
with Ref. iTPTl ). For example, the energy per particle of the 
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FIG. 3: (color online) Spectrum as a function scattering length for 
3 < iV < 6. The thin solid curve corresponds to the trimer state 
energy, the short dashed curves correspond to the four-body states, 
the long-dashed curves correspond to the five-body states, and the 
thick solid curves correspond to the six-body states. Inset: States 
formed below the trimer threshold. 



lowest trimer state of V a at unitarity is ~ 0.0AE sr , while the 
energy per particle of the lowest six-body state at unitarity is 
- 0.6E sr . The latter result implies that Eg /Eg ~ 30, al- 
most a factor of two larger than the universal predictions. The 
introduction of a repulsive three-body force, as proposed in 
Ref. 1151 . leads to a ratio of Eg /Eg ~ 18, which is signifi- 
cantly closer to the universal predictions. 

The increasing importance of nonuniversal corrections as 
N increases is also reflected in the pair distribution functions 
presented in Fig. |4] As N increases, the lowest cluster states 
become more localized in the nonuniversal region (r < ro) 
and, therefore, become less universal. The single peak struc- 
ture of the pair distribution indicates that these states are ba- 
sically droplets which are mainly described by configurations 
at which all particles are at similar distances. 
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FIG. 4: (color online) Pair distribution functions of the lowest N- 
body cluster states (N = 3 — 6). Inset: Distribution functions for the 
trimer ('3'), the trimer-atom four-body state ('3+1'), and the trimer- 
atom-atom five-body state ('3+1+1'). 

The excited states (with energies £4,2, ^5,2, and £^,2) are 
much closer to the lowest fragmentation threshold and can be 
qualitatively described as cluster-atom "dimers" with one par- 
ticle loosely bound to an A— 1 cluster state. This structure can 
also be identified in the pair correlations that coincide at small 
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r with the N — 1 cluster pair correlation but it has a longer tail 
that describes the cluster-atom correlation (cf. inset in Fig. [4}. 

We also identify another five-body resonance, presented in 
Fig. [3] as -E5.3, that is a state energetically below the lowest 
trimer-fragmentation threshold that can decay into the lowest 
tetramer-atom channel. The energy and the pair-distribution 
function show that the state is qualitatively described as a 
trimer weakly bound to two atoms forming a trimer-atom- 
atom state. The energies at unitarity of such states are 1.3E% 
(V a ) or 1.2E% (V t ), i.e., slightly below the trimer and the ex- 
cited tetramer energies, suggesting that most of the contribu- 
tion of the energy comes from the bonding of the trimer sub- 
cluster. Furthermore, the pair-distribution function shows two 
clear peaks that can be identified as coming from atom-atom 
correlations inside the trimer subcluster and atom-atom corre- 
lations between an atom inside the trimer and an atom outside 
the trimer. At small r, the five-body P5 agrees well with the 
trimer P3, suggesting that the small r part of the P5 comes ex- 
clusively from the trimer subcluster contribution. The relative 
strength between the pair-correlation peaks is consistent with 
this interpretation. 

These states can be experimentally observed in ultracold 
gases through the analysis of A-body recombination pro- 
cesses [22 1 . Current experiments with Cs and Li have iden- 
tify three- and four-body resonances through the observation 
of losses at the predicted resonance positions (3}j6). At low 
temperatures, the A-body resonant enhancement of losses oc- 
curs at the critical interaction strengths at which an A-body 
cluster becomes resonant with the free particle-scattering con- 
tinuum. If the A-body clusters behave universally, the posi- 
tions of the resonances are given by critical scattering lengths 
a* N , which are only controlled by the three-body parameter; 
the ratio between any two a* N is a universal number. We ex- 
plore the critical scattering-length ratios for the lowest cluster 
states in a range of model interactions and find that 0.45 < 
al/a* < 0.47 and 0.63 < a^/a* < 0.67. The description of 
the six-body state for strong repulsive core potentials is chal- 
lenging and does not allow an accurate estimation of a%. In 
the range of potentials where convergence is achieved we find 
that al/a* b ~ 0.73 - 0.74. The four-body predictions are 
relatively close to the universal prediction of a\/al « 0.43 
and its deviations are comparable to those observed in experi- 
ments (which also analyze the first Efimov family). The five- 
and six-body scattering-length ratios are in the same ballpark 
of Ref. ED predictions (a* 5 /al w 0.69 and a* 6 /a* 5 « 0.78) 
and of Ref. [B) (a* 5 /a* 4 « 0.6 and a* 6 /a* 5 w 0.7) based on 
different model interactions. From the analysis of the five and 
six-body resonances we predict universal scattering-length ra- 
tio of al/a\ w 0.66 and a* 6 /a* 5 0.78. 

In conclusion, we have demonstrated the existence of uni- 
versal five- and six-body resonances. In experiments, these 
five- and six-body states should manifest as a loss peak at the 
critical scattering length that should occur at approximately 
ag 6 ~ 270ao and at ag b ~ 210cto for the Cs experiments at 
Innsbruck f3] IU. Our results strongly support the hypothe- 
sis of a universal regime in which bosonic systems are only 



controlled by two- and three-body physics. We extend the 
analysis to the lowest Efimov family and we identify a rich 
structure of bound states and resonances in bosonic systems 
with large scattering lengths. This structure can be interpreted 
as emerging from a universal behavior that is strongly affected 
by nonuniversal corrections. Therefore, much of this structure 
is expected to persist in the universal regime. In particular, 
we expect the persistence of universal trimer-atom-atom res- 
onant states. A recent study ifTTI predicts a large and positive 
atom-trimer scattering length, which implies that, even in the 
universal regime, the trimer-atom-atom system fulfils Efimov 
conditions for the formation of weakly bound trimers. The 
emergence of universal picture for bosons suggests a reinter- 
pretation of previous studies on strongly interacting systems. 
For example, one can speculate that the similarities found in 
the formation of 4 He and Tritium clusters E4l have their root 
in an underlying universal behavior that has Efimov physics at 
its root, but is modified by nonuniversal corrections. 
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